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A  PNS  code's  capability  to  predict  the  onset 
of  laminar  and  turbulent  boundary  layer 
separation,  forces,  moments,  and  pressure 
distributions  was  Investigated.  Predictions  of 
incipient  separation  show  good  agreement  with 
experimental  results.  Force  and  moment 
calculations  were  found  to  agree  favorably  with 
experimental  data  for  a  wide  range  of  geometries 
for  Mach  numbers  from  2  to  U  with  turbulent 
boundary  layers.  Limited  comparison  of  pressure 
distribution  data  also  showed  good  agreement. 

Skin  friction  and  Stanton  number  calculations  are 
presented  to  indicate  the  code's  capability  to 
calculate  this  Information  and  hence,  provide  a 
complete  data  base  for  attached  flow  conditions. 
The  sensitivity  of  these  calculations  to 
smoothing  and  marching  parameter  variations  was 
examined  and  found  to  be  acceptable  for  most 
design  work.  The  code  was  found  to  be  suitable 
to  complement  or  reduce  wind  tunnel  experiment¬ 
ation  for  the  design  of  this  class  of  configura¬ 
tion  with  regard  to  accuracy,  cost  and  usability. 

Nomenclature 

A  =  reference  area 

Cf  c  skin  friction  coefficient, 

=  pitching  moment  coefficient  about 
model  tip,  (moment  )/1/2  P.u^Ad 

:  normal  force  coefficient, 

( normal  force ) !  1  /2o„„u^A 

Cq  =  initial  pitching  moment  curve  slope, 

‘  ,  per  deg . 

3a 

Ik  =  initial  normal  force  curve  slope, 

.  per  deg. 

•a 

;  pressure  coefficient,  (P-Po,)/l 
-p  :  specific  heal  at  constant  pressure 

^  :  cylinder  diameter 

-  length  of  cone-cylinder  combination. 


-  total  body  length 
-■  Length  of  cone 


•Aerospace  Engineer,  Aeromechanics  Division, 
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=  length  of  cylinder 
=  length  of  flare 
=  freestream  Mach  number 
=  pressure  1 
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:  w-all  heating  rate 


=  Stanton  number, 

0w/P<A>„Op(Tw  -  Tt) 

=  total  temperature 

=  wall  temperature 

:  axial  component  of  freestream  velocity 


:  axial  location 


z  angle  of  attack 
=  cone  half  angle 
=  flare  angle 
:  freestream  viscosity 
=  freestream  density 


=  axial  component  of  wall  shear 
=  circumferential  angle 


I.  Introduction 

There  are  an  abundance  of  technical  papers  in 
the  literature  concerning  the  application  of  PNS 
codes  to  various  configurations;  see,  for  example, 
References  1-7.  A  summary  of  PNS  work  is  given  In 
Reference  8.  Some  of  these  papers  address  the 
flow  about  flaps  and  flares^-S,  Others  address 
the  flow  about  boattalls5“7 ,  However,  little 
effort  has  been  made  to  assess  the  accuracy,  cost 
effectiveness,  and  usability  of  a  PNS  code  for  a 
complete  set  of  data  that  might  be  used  to  design 
a  relatively  simple  aerodynamic  shape,  such  as  a 
cone-cy Under-flared  projectile.  Further,  there 
has  been  little  effort  to  establish  the  limits  of 
a  PNS  code  in  predicting  Incipient  (onset  of) 
separation  with  regard  to  flap  or  flare  design. 

The  parabolic  nature  of  the  PNS  (marching)  codes 
causes  them  to  diverge  or  "blow  up"  when  flow 
velocities  opposite  to  the  marching  direction  are 
encountered3 .  However,  it  has  not  been 
established  that  tht»jd|^rgenoe^orre^oi^s  to 


:e^orresponds  to 

■  i  10  0  -F  9 


Incipient  separation  determined  from 
experimental  data. 


The  pxjrpose  of  this  paper  is  twofold.  The 
first  purpose  is  to  assess  the  ability  of  a  PNS 
code  to  accurately  detect  the  onset  of 
separation.  The  second  purpose  is  to  provide  a 
set  of  calculations  for  cone-cylinder-flared 
bodies  of  revolution,  compare  them  to 
experimental  results,  and  make  an  assessment  of 
the  suitability  of  PNS  calculations  to  supplant 
or  replace  wind  tunnel  testing  including  accuracy 
and  coat  effectiveness.  In  the  next  two  sections 
a  brief  description  of  the  PNS  code  and  the 
configurations  and  flow  conditions  considered  for 
both  the  separation  and  aerodynamic  studies  will 
be  presented.  In  the  third  section,  typical 
results  will  be  examined  and  discussed  for  each 
study,  and  finally,  major  conclusions  will  be 
given . 

II.  PNS  Code  Description 

The  code  chosen  for  this  study  was  developed 
to  solve  the  steady-state,  ideal  gas,  high 
Reynolds  number,  thin  layer,  parabolized  Navler- 
Stokes  equations^"' 1 .  This  marching  code  is 
based  on  a  non-iterative  technique  formulated  by 
Schlff  and  Steger  which  uses  the  Beam  and  Warming 
implicit  finite-difference  relations  and  the 
alternating  direction  implicit  algebraic  equation 

solver^2 ,13 . 

The  initial  code  was  modified  by  Chaussee  et 
al.,2>9-1l  to  solve  for  the  flow  field  around 
configurations  with  blunt  tips  by  using  In-flow 
data  starting  planes  generated  from  a  full  Navier- 
Stokes  code.  The  present  form  of  the  code 
expands  the  work  started  by  Chaussee  et  al.,  to 
include  refinements  in  the  prediction  procedure 
involving  grid  generation,  boundary  conditions, 
and  increased  numerical  stabllity^^.  When 
flows  with  turbulent  boundary  layers  are 
considered,  the  Baldwin-Lomax  algebraic  two  layer 
turbulence  model  is  employed  to  calculate  the 
turbulent  eddy  vlscosity^^. 

All  PNS  flow  field  solutions  were  obtained  by 
marching  from  an  externally  supplied  starting 
solution  one  half  cylinder  diameter  aft  of  the 
conical  nose  vertex.  An  algebraic  grid  measuring 
19  by  45  points  in  the  circumferential  and  radial 
directions,  respectively,  was  employed.  Thirty 
five  of  the  45  radial  grid  points  were  clustered 
near  the  wall  to  resolve  the  boundary  layer.  An 
axial  marching  step  size  of  .025  cylinder 
diameters  was  employed  up  to  a  plane  one  cylinder 
diameter  aft  of  the  cone-cylinder  Juncture.  This 
small  step  size  was  necessary  to  capture  the 
strong  gradients  in  the  flow  field  that  occur  in 
the  shoulder  region.  The  step  size  was  then 
relaxed  to  0.1  diameters  while  marching  along  the 
cylinder  to  a  plane  0.6  diameters  upstream  of  the 
cylinder-flare  Juncture.  The  step  size  was  then 
reduced  to  .025  diameters  to  accommodate  the 
accompanying  rapid  pressure  rise.  This  step  size 
was  approximately  one  fourth  the  boundary  layer 
thickness  in  this  region.  Previous  work  done  by 
Viegas  and  Horstman  indicates  that  a  step  size  on 
the  order  of  half  a  boundary  layer  thickness  is 
sufficient  to  resolve  two-dimensional  shock' 
boundary  layer  lnteraotlons19 .  Nominal  values 
of  the  code's  smoothing  and  damping  parameters 


were  0.2  with  the  exception  of  the  implicit 
smoothing  parameter  which  was  0.40.  The  nominal 
value  of  the  grid  spacing  parameter  at  the  wall 
was  .0001 .  Any  exceptions  to  these  parameters 
will  be  noted  in  the  text. 

III.  Configurations  and  Plow  Conditions 

To  meet  this  paper's  first  objective,  the 
separation  study,  a  10°  cone-cyllnder-fleu'e 
(Pig.  1)  was  examined  at  a  freestream  Mach  number 
of  5  to  determine  the  capability  of  the  PNS  code 
to  detect  the  onset  of  flow  separation  induced  by 
shock/boundary  layer  Interaction  at  the  cylinder- 
flare  Junction.  The  flare  angle  and  Reynolds 
number  were  varied  to  examine  the  relationship 
among  flow  and  geometry  parameters  regarding  the 
onset  of  separation.  Incipient  separation  was 
Investigated  under  both  laminar  and  turbulent 
conditions.  The  laminar  unit  Reynolds  number  was 
varied  from  1.5  x  105/in.  to  4.0  x  105/in. 

The  turbulent  range  was  taken  from  7.0 
X  10^/ln.  to  1 .5  X  10^/in.  All  cases  were 
run  for  a  wall  temperature  of  520°  R  and  a 
IVeestream  temperature  of  432°  R. 

To  address  the  second  objective,  the 
aerodynamic  study,  several  cone-cylinder-flared 
configurations  (Pig.  1)  have  been  examined. 

These  configurations  were  identical  to  those 
tested  in  a  wind  tunnel  study  in  the  mid 
1950s^^“^®.  PNS  results  are  presented  for  Mach 
numbers  of  2.18,  2.6l,  and  4.04  at  unit  Reynolds 
numbers  of  5.0  x  lO^/ln.,  6.0  x  lO^/in.,  and 
4.0  X  105/ln.,  respectively.  Pully  turbulent 
flow  was  assumed  for  all  conditions.  All  cases 
were  run  for  a  freestream  total  temperature  of 
542°  R  and  a  wall  temperature  of  527°  R.  This 
wall  temperature  is  very  near  adiabatic 
conditions.  An  adiabatic  wall  condition  was 
tried,  but  the~code  failed  to  march~^  to  the. 20° 
flare .  however,  the  code  successfuTTy  marched 
over  a  5°  flare  under  these  conditions. 

IV .  Results  and  Discussion 
A.  Separation  Study 

The  separation  study  with  its  objective  of 
determining  the  ability  of  the  PNS  code  to  detect 
incipient  separation  was  accomplished  for  a  given 
Reynolds  number  by  gradually  increasing  the  flare 
angle  until  the  code  indicated  separation  had 
occurred,  and  then  comparing  the  results  with 
experimental  data.  These  results  are  compared  in 
Pig.  2  with  experimental  two-dimensional  ramp 
data  on  incipient  separation  correlated  by 
Needham  and  Stollery‘® .  To  first  order  the 
flare  appears  as  a  two-dimensional  ramp  to  the 
PNS  code.  The  PNS  code  provides  reasonable 
estimates  for  the  condition  of  incipient 
separation.  Por  the  case  of  laminar  flow, 
separation  is  indicated  at  slightly  higher  flare 
angles  than  occurred  for  the  experimental  data  on 
two-dimensional  ramps.  The  trend  with  Reynolds 
number  closely  follows  the  experimental  data. 

Por  turbulent  flow,  separation  occurs  at  a 
slightly  lower  flare  angle  than  the  two- 
dimensional  data  would  Indicate.  The  affect 
on  indicated  separation  caused  by  variations  in 
the  code's  smoothing  parameters  was  examined  and 
is  summarized  in  Table  1 .  The  results  indicate 
that  a  decrease  in  the  incremental  smoothing 


parameters  causes  the  code  to  register  axial 
separation  at  lower  flare  angles.  This  Is 
consistent  with  the  purpose  of  these  parameters 
to  augment  the  axial  velocity  component  so  that 
Isolated  small  negative  axial  velocities  will  not 
arbitrarily  preclude  further  marching. 

Similarly,  results  indicate  that  a  decrease  In 
the  explicit  and  Implicit  artificial  viscosity 
smoothing  parameters  causes  the  code  to  indicate 
separation  at  higher  flare  angles.  It  Is  not 
completely  understood  why  this  latter  affect  is 
present.  In  conclusion,  using  the  nominal  values 
of  the  smoothing  parameters,  reasonable  estimates 
of  Incipient  separation  can  be  obtained. 

B .  Aerodynamic  Study 

Over  fifty  different  cases  were  calculated  to 
verify  the  code's  suitability  to  predict  flow 
field  and  aerodynamic  force  and  moment  variations 
due  to  changing  geometric  parameters  and 
freestream  Mach  numbers.  This  number  of  cases 
was  sufficient  to  obtain  an  estimate  for  the  coat 
of  computing  a  complete  data  matrix  as  In  the 
wind  tunnel  testa^®”'®  and  to  obtain  an 
assessment  of  the  utility  of  the  code. 

Typical  force  and  moment  results  are  shown  In 
Figs.  3-10.  Variations  of  normal  force 
coefficient  and  pitching  moment  coefficient  with 
angle  of  attack  are  compared  with  the 
experimental  data  In  Pigs.  3-**.  Variation  of 
Initial  normal  force  curve  slope  and  Initial 
pitching  moment  curve  slope  with  flare  angle  and 
flare  length  are  compared  with  the  experiment  In 
Figs.  5  and  6.  Variations  of  normal  force  curve 
slope  and  pitching  moment  curve  slope  with 
cylinder  length  and  Mach  number,  as  well  as 
comparison  with  experimental  data,  are  presented 
In  Figs.  7-10.  A  typical  axial  distribution  of 
the  pressure  coefficient  for  both  the  windward 
and  leeward  rays  at  2°  angle  of  attack  Is 
presented  and  compared  with  experimental  data  In 
Fig.  11.  Circumferential  distributions  at  three 
axial  locations  are  given  In  Fig.  12.  Similar 
axial  and  circumferential  distributions  for  skin 
friction  coefficient  and  Stanton  number  for  the 
same  conditions  are  presented  In  Figs.  13-16; 
however,  experimental  data  were  not  available  for 
comparison . 

Figures  3-10  demonstrate  basic  force  and 
moment  results  for  parametric  variations  normally 
encountered  In  a  design  process.  All  of  the 
results  show  good  agreement  with  experimental 
data.  The  absolute  values  of  the  calculations 
are  slightly  higher  than  the  corresponding 
experimental  data  (as  much  as  ten  to  fifteen 
percent  for  large  flare  angles  and  angles  of 
attack).  However,  the  accuracy  achieved  Is 
sufficient  for  most  design  work.  As  one  would 
expect,  larger  flares  produce  larger  normal  force 
curve  slopes  and  pitching  moment  curve  slopes. 
This  Is  evident  In  Figs.  5  and  6  for  both  flare 
length  and  flare  angle.  Generally,  the  affect  of 
a  longer  cylinder  length  Is  to  decrease  CNq 
(Fig.  7).  This  Is  a  result  of  so  called  "lift 
carry  over'ZI  uhich  refers  to  the  significant 
circumferential  pressure  difference  around  short 
cylinders  at  angle  of  attack.  This  pressure 
difference  Is  converted  Into  additional  lift  on 
the  flare  above  that  which  the  flare  would 
generate  If  the  pressure  was  more  uniform  around 


the  cylinder,  which  Is  the  ease  for  the  longer 
cylinders  (  ^cy^®  -  'i-OI*  Hence,  the  curves  In 
Fig.  7  are  flattening  out  with  Increased 
The  affect  of  cylinder  length  on  Is  much 
more  dramatic  and  direct.  Although  there  Is  a 
alight  decrease  In  with  cylinder  length  as 

Just  discussed,  the  primary  affect  on  Is  to 

Increase  the  moment  arm.  Hence,  C_  Increases 
almost  linearly  with  ^y/d  (Fig.  ST?  In 
Figs.  9  and  10  the  Mach  number  trends  are  the 
same  as  one  would  obtain  on  cones  (Ref.  22,  Chart 
8).  Cnq  and  decrease  with  Mach  number 

for  3f  =  5°.  However,  both  coefficients 
Increase  with  Mach  number  for  6f  =  20'^. 

Drag  coefficients  were  not  measured  In 
References  16-18;  however,  a  comparison  with 
experimental  data  at  H„  =  3.0  on  a  20°  flare®3 
shows  good  agreement: 


Cd  ,  2.22 

“exper Iment 


,  M„  =  3 .0 


fig  =  22.5  ,  <cy/d  =  “I.  if/a  -  2.‘* 

Cn  .  =  2.19  ,  M„  =  2.81 

^computation 

This  Is  not  surprising  since  the  total  drag  Is 
dominated  by  wave  drag  (base  drag  Is  not 
included),  and  the  pressure  distribution  shown  In 
Fig.  11,  to  be  discussed  next,  agrees  quite  well 
with  the  experiment. 

Figure  1 1  demonstrates  the  generally  good 
agreement  between  the  computed  and  experimental 
axial  pressure  coefficients,  with  the 
experimental  data  being  slightly  higher  than  the 
computations.  Fig.  12  shows  the  circumferential 
distribution  of  pressure  coefficients.  The 
distributions  immediately  upstream  of  the  flare 
(x/l  s  0.50)  and  at  the  flare's  base  (x/1  =  1.0) 
show  the  expected  steady  decrease  In  pressure 
from  windward  to  leeward,  although  the  variation 
Just  before  the  flare  Is  very  small.  The 
distribution  Just  after  the  cylinder-flare 
Junction  (  x/1  =  0.5*1)  begins  with  this  trend,  but 
the  pressure  starts  to  rise  again  as  It 
approaches  the  leeward  side  ( ♦  =  70°)  and  then 
stabilizes.  This  suggests  the  possibility  that 
circumferential  separation  Is  beginning  to 
occur.  Circumferential  separation.  If  present. 

Is  not  very  strong  since  there  Is  no  evidence  of 
It  further  back  on  the  flare.  (Though  not 
Included  in  the  results,  a  long  cone -cylinder 
configuration  was  examined  at  these  flow 
conditions  and  showed  strong  evidence  of 
circumferential  separation.)  It  should  be  noted 
that  this  20°  flare  angle  Is  very  near  incipient 
separation,  as  can  be  deduced  from  Fig.  2.  This 
Is  supported  by  the  fact  that  the  experimental 
pressure  distribution  (Fig.  11)  also  shows  a 
slight  pressure  rise  starting  at  approximately 
x/ i  -  0.*l8.  The  computed  results  for  skin 
friction  and  Stanton  number  for  this  case  also 
Indicate  that  the  boundary  layer  Is  close  to 
Incipient  axial  separation. 

The  axial  and  circumferential  distributions 
of  the  skin  friction  coefficient  are  presented  In 
Figs.  13  and  11,  respectively.  An  examination  of 
Fig.  13  reveals  that  skin  friction  varies 
smoothly  with  axial  position  except  In  the 
vicinity  Just  downstream  of  the  flare  Junction. 


iTi^irv  L-W  ‘ 


I  ** 


The  flow  expansion  at  the  cone/cy Under  shoulder 
Increases  the  boundary  layer  thickness  resulting 
In  a  rapid  decrease  in  skin  rrlctlon  coefficient 
which  then  steadies  to  a  slow  decline  as  one 
expects  of  a  turbulent  boundary  layer.  The 
sudden  increase  in  the  skin  friction  coefficient 
at  the  flare  Junction  is  physically  unreal  and  is 
a  result  of  the  parabolic  nature  of  the 
equations.  In  the  real  flow  the  pressure  rise  is 
felt  forward  of  the  flare  Junction,  allowing  the 
flow  to  gradually  decelerate  reducing  the  skin 
friction  ahead  of  the  flare  with  a  further  rapid 
decrease  on  the  first  position  of  the  flare.  The 
existence  of  the  negative  value  of  akin  friction 
can  be  misleading  since  the  code's  smoothing 
parameters  produce  the  necessary  stabilization 
under  these  conditions.  However,  it  does 
indicate  that  the  20°  flare  la  very  close  to 
axial  separation.  The  flow  finally  stabilizes, 
and  the  akin  friction  coefficient  on  the  flare 
resumes  its  steady  decline.  As  expected,  the 
skin  friction  is  higher  here  than  on  the  cylinder 
because  of  the  thinner  boundary  layer  on  the 
flare.  The  circumferential  distributions  at 
x/1  =  0.50  and  x/1  =  1.0  vary  smoothly, 
reflecting  the  thinner  boundary  layer  on  the 
windward  aide  (higher  Cf)  and  thicker  boundary 
layer  on  the  leeward  side  (lower  Cf ) ;  see  Fig. 

1*1.  The  distribution  at  x/2.  =  O.si  la  not  as 
uniform  and  reflects  the  fact  that  there  la  some 
circumferential  separation  occurring  Immediately 
downstream  of  the  flare  Junction.  Details  of 
this  curve  and  similar  curves  in  the  region 
between  x/1  =  0.52  and  0.60  should  not  be  taken 
as  absolute  but  do  indicate  boundary  layer  trends 
and  the  nearness  of  axial  separation. 

Figures  15  and  16  respectively  show  the  axial 
and  circumferential  Stanton  number.  The  computed 
values  result  from  the  wall  temperature  being  set 
slightly  higher  than  the  wall's  recovery 
temperature  under  adiabatic  conditions 
(Ty  =  527°  R) .  Hence,  In  general  the  heat 
transfer  is  from  the  wall  to  the  flow  which 
results  In  a  negative  Stanton  number  over  much  of 
the  body  (Fig.  15).  Exceptions  occur  downstream 
of  the  cone-cylinder  shoulder  and  downstream  of 
the  flare  Junction.  At  the  shoulder,  the 
boundary  layer  is  cooled  by  the  expansion  and 
greatly  thickened.  The  cooling  would  tend  to 
make  the  heat  transfer  more  negative ,  and  the 
thickening  would  reduce  the  magnitude.  Neither 
would  change  the  sign  as  Indicated.  Therefore, 
the  change  In  sign  Is  not  fully  understood  and 
may  be  a  result  of  the  parabolic  nature  of  the 
equations.  Downstream  of  the  flare  Junction  the 
compression  heats  the  flow  Increasing  the  air 
temperature ,  and  the  boundary  layer  becomes 
thinner  due  to  the  pressure  rise.  The 
compression  Is  sufficient  to  change  the  sign  of 
the  heat  transfer.  The  thinning  of  the  boundary 
layer  would  increase  the  magnitude.  Neither, 
however,  are  sufficient  to  explain  the  very  high 
Stanton  number  Immediately  after  the  flare 
Junction.  Again,  this  spike  may  be  a  result  of 
the  parabolic  nature  of  the  equations.  Further 
downstream  of  the  Junction,  the  flow  rapidly 
adjusts  and  resumes  its  negative  Stanton  number 
level.  The  circumferential  variation  tends  to  be 
smooth  at  the  * '  <  =  0.50  and  1.0  locations,  being 
slightly  more  negative  on  the  windward  side  where 
the  boundary  layer  Is  thinner  (Fig.  16). 

However,  at  x  '  r  0.5U  a  change  In  the  -Stanton 


number  near  t  -70°  is  evident.  Again,  this  is 
associated  with  a  possible  circumferential 
separation. 

A  very  brief  study  was  conducted  to  determine 
the  sensitivity  of  results  to  code  Inputs  such  as 
smoothing  parameters,  marching  step  size,  and 
grid  spacing  in  the  boundary  layer.  The  affect 
of  smoothing  parameters  on  axial  separation  has 
been  discussed  previously.  The  force  and  moment 
results  were  only  weakly  affected  by  these 
parameters,  on  the  order  of  a  few  percent  or 
less.  All  the  distributions  were  likewise  only 
slightly  affected  except  in  the  region  of 
shook/boundary  layer  interaction  near  the 
Junction  of  the  flare,  partlcualarly  the  20° 
flare  angle  case  where  the  flow  is  near 
separation.  In  this  area  significant  affects 
were  present.  The  pressure  coefficient  varied  by 
thirty  to  forty  percent  near  the  maximum  Cp, 
and  the  oscillations  that  occurred  did  not  dampen 
out  immediately.  The  greatest  sensitivity  due  to 
the  variation  of  these  parameters  occurred  with 
the  skin  friction  coefficient  in  this  same 
region.  This  is  not  as  critical  in  one  sense, 
because  skin  friction  is  not  an  Important 
contribution  to  drag.  However,  the  skin  friction 
curves  may  be  the  best  indicator  of  how  close  to 
separation  the  flow  is,  rather  than  relying  on 
the  PNS  code  falling  as  a  criteria.  To 
summarize,  the  sensitivity  study,  although  brief, 
indicated  that  basic  design  data,  such  as  forces 
and  moments,  are  insensitive  to  the  above 
parameters  as  are  the  distributions  of  pressure, 
skin  Criotlon,  and  Stanton  number,  except  in  the 
region  of  a  rapid  geometry  change.  However,  in  a 
region  of  rapid  geometry  change  such  as  the  flare 
Junction,  results  are  sensitive  to  these 
parameters,  and  the  user  should  conduct  his  own 
sensitivity  studies. 

C,  Evaluation  of  Usability 

It  is  not  sufficient  that  the  PNS  code 
provide  accurate  results.  The  code  must  do  so  in 
a  cost  effective  fashion  as  well  as  provide 
results  in  a  timely  and  flexible  manner  to  impact 
the  engineering  design  of  an  airframe.  In 
discussing  the  cost  and  usability  of  the  code, 
two  sets  of  factors  must  be  remembered.  First, 
the  results  were  generated  on  a  medium  size 
computer,  a  CDC  Cyber  176.  Second,  the  users,  in 
this  case  the  authors,  are  not  computational 
experts  but  users  of  codes  for  aerodynamic 
research  and  analysis  along  with  other  available 
tools,  both  experimental  and  analytic. 

It  is  estimated  that  a  complete  computational 
data  matrix  comparable  to  that  of  the  wind  tunnel 
test  of  Ref.  16-18  could  be  obtained  at  one  sixth 
the  cost  of  a  wind  tunnel  test  today.  This 
computational  estimate  includes  a  small  factor  to 
cover  geometry  initiation  and  errors  in  data 
entry.  The  computational  cost  also  reflects 
maximum  use  of  restarts  using  stored  data  planes 
generated  from  previous  calculations.  This 
factor  of  six  is  considered  to  be  a  conservative 
estimate,  because  the  low  computer  cost, 
flexibility,  and  short  response  time  (on  the 
order  of  a  few  hours)  for  the  computation  would 
allow  a  more  efficient  approach  to  the  design 
process.  With  this  flexibility  and  short 
response  time  only  a  skeleton  matrix  would  be 
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coBputed  to  determine  basic  trends  such  as  have 
been  Illustrated  here.  A  preliminary  design 
could  be  developed  from  these  trends,  and  from 
then  on,  only  limited  computations  would  be 
performed  as  the  design  matured.  Because  of 
scheduling  constraints  this  approach  la  nearly 
Impossible  with  a  wind  tunnel.  One  final  factor 
that  needs  to  be  considered  In  the  coat  analysis 
Is  that  the  calculations  provide  not  only  the 
forces  and  moments,  but  the  entire  flow  field 
Including  skin  friction,  heat  transfer  and 
pressure  distributions  at  no  additional  coat. 
Normally,  separate  tunnel  entries  would  have  to 
be  made  for  each,  or.  In  many  cases  the  designer 
might  have  to  do  without  skin  friction,  heat 
transfer,  or  pressure  distribution  results.  This 
additional  Information  can  lead  to  a  more 
complete  design  and  one  In  which  there  Is  more 
confidence  because  of  that  completeness.  When 
all  of  these  factors  are  taken  into  account  the 
cost  of  design  by  computation  of  these  simple 
bodies  with  attached  flow  la  an  order  of 
magnitude  leas  costly  than  for  wind  tunnel 
testing.  In  addition,  as  Is  noted  in  Reference 
2U ,  the  two  factors  of  flexibility  and  Increased 
confidence  are  of  more  importance  than  cost  In 
many  cases.  These  factors  Impact  other  areas  of 
the  design  and  hence  can  reduce  the  overall 
design  cycle  time  and  thus  reduce  overall  coat. 

V.  Concluding  Remarks 
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9.  Chaussee,  D.S.,  Kutler,  P.,  and  Pulliam,  T. 
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Program:  Part  1 ,  Viscous  Blunt  Body  Program 
(Interim  Report),"  AFWAL-TM-81-63-FIMG,  Mar  8l. 

10.  Chaussee,  D.S.  and  Steger,  J.L., 
"Three-Dimensional  Viscous  Flow  Field  Program: 
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Based  on  calculations  with  a  PNS  code  of  a 
large  number  of  cone-cylinder -flared 
configurations  at  several  Mach  numbers  and 
Reynolds  numbers  under  attached  flow  conditions 
the  following  conclusions  can  be  drawn. 

The  PNS  code  examined  provided  flow  field  and 
aerodynamic  force  and  moment  data  of  sufficient 
accuracy  for  general  design  of  this  class  of 
configurations.  The  PNS  code  provided  reasonable 
estimates  for  the  onset  of  separation  on  flared 
bodies  and  hence,  can  be  used  as  a  tool  to 
establish  the  design  limits  for  attached  flow. 

The  PNS  calculations  are  a  very  cost  effective 
approach  for  the  design  of  cone-cy Under-flared 
configurations  with  attached  flow. 
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Fig.  2  Correlation  of  Incipient  Separation. 
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Fig.  3  Norsal  Force  Coefficient  va  Angle  of 
Attack. 


Fig.  1  Tone-Cy llnder-Klare  Geometry 
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Fig.  S  Initial  Normal  Force  Curve  Slop*  v# 
Flare  Length. 


rig.  7  Initial  Normal  Force  Curve  Slope 
vs  Cylinder  Length. 


V/  ■ •  n.s-.Of  -  20- 


ri.s-  Ifid  ■ ‘i.h.  i,y/d 

o  OATAtRCfS  ItrlS) 

a  COAi^UTfP  POIM11 


O  ci^catNCNUi.  DAfAUfr  lb. 
Q  COn^UTCD  POINTS 


Fig.  8  Inlttal  Pitching  Moment  Curve  Slope 
vs  Cylinder  Length. 
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Fig.  10  Initial  Pitching  Moaent  Slope  va 
Mach  Number. 
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Fig.  11  Pressure  Coefficient  vs  Axial 
Position. 
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g  Initial  Normal  Force  Curve  Slope  va 
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Fig.  12  Pressure  Coefficient  vs 
Circumferential  Angle. 


